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A fundamental understanding of phonon transport in stanene is crucial to predict 
the thermal performance in potential stanene-based devices. By combining hrst- 
principle calculation and phonon Boltzmann transport equation, we obtain the lattice 
thermal conductivity of stanene. A much lower thermal conductivity (11.6 W/mK) 
is observed in stanene, which indicates higher thermoelectric efficiency over other 
2D materials. The contributions of acoustic and optical phonons to the lattice ther¬ 
mal conductivity are evaluated. Detailed analysis of phase space for three-phonon 
processes shows that phonon scattering channels LA-l-LA/TA/ZAoTA/ZA are re¬ 
stricted, leading to the dominant contributions of high-group-velocity LA phonons 
to the thermal conductivity. The size dependence of thermal conductivity is investi¬ 
gated as well for the purpose of the design of thermoelectric nanostructures. 

I. INTRODUCTION 

Research in 2D materials has experienced an explosion of interest due to their potential 
for integration into next-generation electronic and energy conversion deviced^. Recently, a 
new 2D material, stanene, has been realized on the Bi 2 Te 3 (lll) substrate by molecular beam 
epitax}^, and attracted tremendous interest due to its extraordinary properties. Stanene, 
a monolayer of tin him with buckled honeycomb structure like silicene, is predicted to be 
an example of a topological insulator, which can support a large-gap 2D quantum spin Hall 
(QSH) statd^. It has been reported that a nontrivial bulk gap of ~0.1 eV opens when the 
spin-orbital coupling (SOC) is present in stanene, and the edge states are gapless with bands 
dispersing inside the bulk gap and helical with spin-momentum locking, which can be used as 
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dissipationless conducting “wires” for electronic circuitd^. The free standing stanene with 
~0.1 eV SOC gap is a promising applicant for room temperature and even high temperature 
application. In addition, a recent study has found that the thermoelectric (TE) hgure of 
merit zT is strongly size dependent in stanene, and can be improved to be signihcantly 
larger than 1 by optimizing the geometric size, which provides great opportunities in the 
next-generation TE device^. 

All these electrical and thermoelectrical applications of stanene are closely related to its 
thermal transport properties. For instance, high TE performance requires the system to 
be a bad conductor for phonons but a good conductor for electrons, which often conflict 
with each othei^^^^. In 2D QSH systems such as stanene, the edge states form perfectly 
conducting channels for electrons, leading to great improvement in the TE hgure of merit 
;^7 lio | i3|i4 ]_ Moreover, zT of stanene can be further improved by optimizing the geometry 
size to maximize the contribution of the gapless edge states, which signihcantly helps to 
optimize the electron transport. However, the phonon transport properties of stanene have 
not been investigated in detail. Furthermore, although the size dependence of zT and the 
Seebeck coefficient in stanene has been studied for TE applications, the size dependence 
of lattice thermal conductivity, which is also crucial when designing TE nanostructures, 
remains uninvestigated. 

Generally both electrons and phonons contribute to the total thermal conductivity. In 
this paper only the lattice thermal conductivity of stanene is under consideration. We cal¬ 
culate the lattice thermal conductivity k of stanene using hrst-principle calculation and an 
iterative solution of the Boltzmann transport equation (BTE) for phononJ^. A much lower 
thermal conductivity (11.6 W/mK) is observed in stanene compared to other 2D materials, 
which makes stanene a promising candidate to realize high thermoelectric performance. We 
also extract the contribution of each phonon mode to investigate the underlying mechanism 
behind the novel thermal transport properties of stanene, and hnd that LA phonons con¬ 
tribute most to the thermal conductivity (>57%) over a wide temperature range. It can 
be attributed to their high phonon group velocities and restricted phase space for three- 
phonon processes. Finally, we obtain the representative mean free path (MFP) of stanene 
that is useful when designing TE nanostructures. The MFP in the small-grain limit is 
also obtained when the nanostructuring induced phonon scattering dominates. The role of 
boundary scattering in stanene nanowires is examined as well. 
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FIG. 1: (a)Top view and (b)side view of stanene. 



FIG. 2: Phonon dispersion and DOS of stanene. Three acoustic and three optical phonon modes 
are indicated with different colors. 

II. RESULTS AND DISCUSSIONS 

A. Crystal structure and phonon dispersion 

Fig. ID shows the optimized structure of stanene. A similar low-buckled conhguration is 
also found in silicene, germanene and blue phosphorend®®®^, which is different from the 
planar geometry of graphend^. The optimized lattice constant and buckling height are 4.67 
A and 0.85 A, respectively, which are in good agreement with the previous workP®. There 
are two tin atoms in one unit cell, corresponding to three acoustic and three optical phonon 
branches. 

Fig. I shows the calculated phonon dispersion. Similar to a-Sn and /3-SrP^, the longi¬ 
tudinal acoustic (LA), transverse acoustic (TA), and z-direction acoustic (ZA) branches of 
stanene are linear near the T point, which is distinctly different from the quadratic ZA- 
branch dispersion in graphend^. This difference between graphene and stanene comes from 
structural difference. In planar graphene, the quadratic ZA-branch dispersion near the T 
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point can be explained as the consequence of the complete decoupling of the Z modes and 
XY modes due to the Dq^ point-group symmetrj^^IM] However, the buckled structure of 
stanene leads to the point group symmetry of which breaks the reflection symmetry 
{z —)■ —z), and therefore the Z modes of stanene couple with the XY mode^, which results 
in the linear dispersion of ZA phonons near the F point. 

Using the phonon spectrum, we calculate the group velocities of acoustic phonons along 
F-M and F-K respectively, as shown in Fig. [^a) and (b). The anisotropy in group velocities 
along the two directions is weaker than that of black phosphorene and SnSd^^^. As shown 
in Fig. [^c), the phonon group velocities indicate weak anisotropy within the entire hrst 
Brillouin zone (BZ). Much smaller sound velocities are found in stanene (1300-3600 m/s 
compared to 5,400-8,800 m/s in silicend^SI, 4200-6,800 m/s in MoS^, 4,000-8,000 m/s in 
blue phosphoren^, and 3,700-6,000 m/s in graphend^. For LA, TA and ZA phonons, 
the sound velocities at long-wavelength limit are about 3548.7 m/s, 1811.4 m/s and 1303.7- 
1536.8 m/s, respectively. 

The Griineisen parameter , which provides information on the anharmonic interactions, 
is given by 


7i(q) = 


ao duj{q) 


( 1 ) 


Wj/q) da 

where Oq is the equilibrium lattice constant, j is phonon branch index, and q is wave vector. 
Fig. I^shows the Griineisen parameters within the hrst BZ. Strong anisotropy is found in the 
Griineisen parameters of TA phonons at high frequencies, while the isotropy in the Griineisen 
parameters of LA and TA phonons at low frequencies indicates that the anharmonicity 
of these phonons is similar from one phonon wave vector to another. The low-frequency 
phonons contribute most to the lattice thermal conductivity of stanene due to their higher 
group velocities and weaker scattering. Therefore, the orientation-independent phonon group 
velocities and Griineisen parameters at low frequencies are crucial to understand the phonon 
transport in stanene. 


B. Lattice thermal conductivity 

The intrinsic lattice thermal conductivity k of naturally occurring stanene is calculated 
using both the iterative and single-mode relaxation time approximation (SMA) solution, as 
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FIG. 3: Phonon group velocities along the (a)r-M and (b)r-K direction, and (c)within the first 
BZ. 

shown in Fig. In contrast to the anisotropy of thermal conductivity observed in black 
phosphorene and which possess the hinge-like structure, the isotropic k is found 

in stanene, which is similar to other 2D hexagonal materials with symmetry such as 
silicene, blue phosphorene, graphene, and monolayer 

The thermal conductivity from the iterative solution at 300 K is 11.6 W/mK and 11.9 
W/mK for naturally occurring and isotopically pure stanene respectively. Naturally occur¬ 
ring Sn consists of 0.97% ^^^Sn, 0.65% ^^^Sn, 0.36% ^^®Sn, 14.7% ^^®Sn, 7.7% ^^'^Sn, 24.3% 
ii^Sn, 8.6% 32.4% ^^oSn, 4.6% i22Sn, and 5.6% 

The measured room-temperature k of 500-nm-thick tin him is 46 ± 4.2 W/mK while for 
the thickness of 100 nm, k is 36 ±2.88 W/mKl^, which are lower than the room-temperature 
K of bulk tin, 64 W/mK. Note that when the thickness of the Sn thin him decreases, phonons 
become more sensitive to boundary scattering®, which further reduces the value of k. Thus 
a smaller calculated k, for stanene is a reasonable expectation. The Debye temperature 0 d 
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FIG. 4: Calculated Griineisen parameters within the first BZ. 
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FIG. 5: Calculated lattice thermal conductivity k of stanene as a function of temperature ranging 
from 100 K to 800 K using the iterative (red) and SMA (blue) solution of the BTE. The inverse 
relationship of k and temperature (k ~ 1/T) is displayed as dot lines. 

can be calculated from the highest frequency of normal mode vibration (Debye frequency) 
I'm = 1-51 THz, 


©D = hum/ks (2) 

where h is the Planck constant, and is the Boltzmann constant. The calculated Debye 
temperature is 72.5 K, which is nearly three times lower than that of /3-Sn of 200 ± 3 KI^SI. 
Since the Debye temperature is a measure of the temperature above which all modes begin 
to be excited and below which modes begin to be frozen oulP^i^, this low ©/) in stanene 
also indicates that the k of stanene is much lower than that of bulk tin. 
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According to the kinetic formula for the thermal conductivity of an assembly of total 
specihc heat C*, average phonon velocity n, and MFP I, i.e. k = l/3Cvl, when the tem¬ 
perature is well above the Debye temperature Qo, fhe C is constant, while the dominant 
three-phonon anharmonic processes at high temperatures lead to I oc yyg the 

calculated k above Qd, which reveals an inverse relation to temperature, i.e. k ~ 1/T. As 
shown in Fig. the htted k (plotted as dot lines) coincides perfectly with the calculated 
K. Such inverse relationship of k and T is also found in other materials such as graphene, 
black phosphorene, PbSe, and Mg 2 (Si,Sn) when T > 0j2iIi2lS21. Furthermore, we compare 
the thermal conductivity calculated from the iterative and SMA solution. The difference 
between the two approaches is very small. It should be noticed that the SMA assumes that 
individual phonons are excited independently, which has no memory of the initial phonon 
distribution, therefore the SMA is inadequate to describe the momentum-conserving charac¬ 
ter of the Normal processes and it works well only if the Umklapp processes dominat^^. Our 
results indicate the dominant role of the Umklapp processes in phonon-phonon interactions 
over a wide temperature range. 

Compared to other 2D hexagonal materials, the k of stanene is smaller than that of silicene 
(26 W/mK)l®l^ monolayer M 0 S 2 (34.5 ± 4 W/mK)l33Iill ^^d blue phosphorene (78 W/ m kR 
while at least two orders of magnitude smaller than that of graphene (3000 W/mK)E2Eil_ 
The lower k of stanene is due to lower Debye temperature (72.5 K compared to 500 K 
for monolayer M 0 S 2 and blue phosphoren d^^*^^ ( and 2,300 K for graphend^. The Debye 
temperature reflects the magnitude of sound velocity. Lower Debye temperature results in 
decreased phonon velocities as mentioned above, and higher phonon scattering rates since 
more phonon modes are active at a given temperatur d^^ * ^^ *^. 

High TE performance requires the system to be a bad conductor for phonons and a good 
conductor for electrons. As a 2D topological insulator with lower lattice thermal conduc¬ 
tivity, stanene can realize much higher thermoelectric efficiency than other 2D materials, 
which makes it a promising candidate for next-generation thermoelectric devices. 

Furthermore, we examine the contributions of different phonon branches to k, of stanene, 
as shown in Fig. The contribution of LA phonons is more than 57% at temperatures 
ranging from 100 K to 800 K, which is larger than 28% in MoSj^, 26% in blue phosphorend^ 
and 8% in graphend^. The contribution of ZA phonons to the k of stanene (13%) is much 
smaller than that of graphene (80%)!^, blue phosphorene (44%)l43, and monolayer M 0 S 2 
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FIG. 6: Percentage contribution of each phonon mode to k as a function of temperature. 

(39%)E^, but larger than silicene (7.5%)^^. It has been reported that the large contribution 
of the ZA mode to the k of graphene is due to a symmetry selection rule in one-atom-thick 
materials, which strongly restricts anharmonic phonon-phonon scattering of the ZA mod^^, 
while the buckled structure of stanene, silicene, and blue phosphorene breaks out the out- 
of-plane symmetry, in which the selection rule does not apply. As a result, the predicted 
contribution of ZA phonons to k in stanene is much smaller than that in graphene at 300 K. 
To provide more physical insight, we investigate the phonon scattering mechanism in detail. 

C. Relaxation time and phase space for three-phonon processes 

Considering the significant difference in the contribution of acoustic phonon modes to 
the K in stanene and other 2D hexagonal materials, it is worthwhile to investigate the 
relaxation time of each phonon mode as a function of frequency, as shown in Fig. In 
Fig.0 a), the three-phonon relaxation time of LA mode is the longest, and the relaxation 
time of the optical phonons is much shorter than the acoustic phonons. We also calculate the 
isotopic relaxation time, as shown in Fig. [^b). Isotopic scattering processes are wavelength- 
dependent. The long-wavelength phonons can transport all the heat with very little isotopic 
scatterin^^, thus we observe the relative long isotopic relaxation time of acoustic phonons 
at low frequencies. The inverse of total relaxation time is a sum of contributions from 
anharmonic three-phonon scattering and isotopic scatterin^^SEHl^ 


i/D(q) = VD(qr" + i/D(qr- 


( 3 ) 
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FIG. 7: The (a)three-phonon, (b)isotopic and (c)total relaxation time of each phonon branch as a 
function of frequency. 

Thus the three-phonon processes with much shorter relaxation time will dominate the total 
relaxation time. The total relaxation time of stanene is shown in Fig. [^c). It is in the same 
order of magnitude as that of graphene, M 0 S 2 , and silicen^^SIlSISo]^ Fig. | 3 c), the 

total relaxation time of LA mode has large values near the acoustic-optical frequency gap, 
while it is still comparable to that of other acoustic phonon modes at low frequencies. To 
investigate the phonon scattering mechanism more deeply, we calculate the allowed phase 
space for three-phonon processes P 3 . 

In the temperature range where three-phonon processes are dominant, the total phase 
space for three-phonon processes P 3 is dehned 
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where hi is a normalization factor, and 

d*’ = E(5) 

j 

and 

±±q'-G)) (6) 

j'j” 

where Zi)j^^(q) corresponds to absorption processes, i.e. ujj{q) + ojji{q') = ujj//{q + q' — G), 
whereas Zi)j~^(q) corresponds to emission processes, i.e. ujj{q) = + ujj"{q — q' — G). 

According to Eq. (|4||^, P 3 contains a large amount of scattering events that satisfy the 
conservation conditions and can be used to assess quantitatively the number of scattering 
channels available for each phonon mode. Less restricted phase space for three-phonon 
processes implies a larger number of available scattering channels. Consequently, there 
is an inverse relationship between P 3 and the intrinsic lattice thermal conductivity of a 
materiapS^. 

The P3 for absorption and emission processes are shown in Fig. [^a) and (b), respectively. 
The phase space for absorption processes of LA phonons is much smaller than P 3 ^^ 
of ZA and TA phonons, while the phase space for emission processes of LA phonons 
is comparable to P^~^ of TA phonons at low frequencies but much larger than that of ZA 
phonons. According to Eq. Q, the calculated total three-phonon phase space P 3 for LA 
phonons is smaller than P 3 of ZA/TA phonons in a large part of the low-frequency range, 
shown in Fig. [^c). The reason for the relatively small P3 of LA phonons is due to the fact 
that the energy of LA phonons are relatively high and subsequently the scattering channels, 
LA-I-LA/TA/ZA-h-TA/ZA, are restricted. As a result, the relatively small P3 alongside 
with relatively large group velocities of LA phonons as mentioned above, determines that 
LA phonons contribute dominantly to k, of stanene. 

As for ZA phonons, while the relatively low energy of ZA phonons severely restricts the 
corresponding phase space for emission processes P 3 \ it also provides a larger number 
of available scattering channels, ZA-l-TA/ZAoLA/TA, which leads to larger of ZA 
phonons, and hence a reduced contribution to k at low frequencies. 
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FIG. 8: Frequency-dependence of three-phonon scattering phase space at 300 K for (a)absorption 
processes, (b)emission processes and (c)all three-phonon processes. 

It is well known that the dominant contribution of ZA phonons to the k of graphene is 
due to a symmetry selection rule in such one-atom-thick materials, which strongly restricts 
anharmonic phonon-phonon scattering of the ZA mod^®'. The selection rule in graphene 
that arises from the reflection symmetry perpendicular to the graphene plane restricts the 
participation of odd number of ZA phonons in three-phonon processes, e.g. ZA-|-ZA<h-ZA, 
ZA-t-LA/TAoLA/TA. As a result, 60% of the phase space of ZA phonons is forbidden by the 
selection rule, which leads to the dominant contribution of the ZA mode to k in graphene. 
However, for the case of buckled stanene considered here, the reflection symmetry is broken, 
which means the selection rule does not apply. Consequently, the relatively large P 3 of ZA 
phonons leads to less contribution from ZA phonons to k of stanene, i.e. 13% in stanene, 
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compared to 80% in graphene at 300 K. 

D. Size dependence of k 

Since the thermoelectric hgure of merit zT can be improved by optimizing the geometry 
size to maximize the contribution of the gapless edge states, it is also essential to investi¬ 
gate the phonon transport in stanene nanostructures and achieve the optimized transport 
properties by reducing thermal conductivity through nanostructuring in the practical design 
of TE nanodevices. We investigate the size dependence of k by calculating the cumulative 
thermal conductivity with respect to the maximum MFP allowed. The cumulative thermal 
conductivity of stanene at different temperatures are plotted in Fig. |^a). The cumulated k 
keeps increasing as MFP increases, until reaching the thermodynamic limit above a length 
Ldiff which represents the longest mean free path of the heat carrier^^^ES^ 'jpg at 150 
K, 300 K, 450 K, 600 K and 750 K are 107.2 pm, 50.9 pm, 35.1 pm, 25.5 pm and 24.2 
pm, respectively. It is found that phonons with MFPs below 73.9 pm at 150 K, 16.9 pm 
at 300 K, 9.5 pm at 450 K, 8.4 pm at 600 K and 6.6 pm at 750 K, respectively, contribute 
around 75% of the total k, which indicates that the design of nanostructures can be utilized 
to reduce k and enhance zT. In order to properly reduce k for thermoelectric applications, 
we can £t the cumulative k to a uniparametric functiorP^ 

< U.) = (7) 

^max 

where i^max is the ultimate cumulative thermal conductivity, Imax is the maximal MFP con¬ 
cerned, and Iq is the parameter to be evaluated. The htted curves at different temperatures 
are plotted in Fig. [^a), reproducing both the slope and the position of the calculated data 
well when l^ax > h- It yields the parameter /q of 5.72 pm, 2.71 pm, 1.77 pm, 1.31 pm and 
0.95 pm for 150 K, 300 K, 450 K, 600 K and 750 K, respectively, which can be interpreted 
as representative of the MFP of relevant heat-carrying phonons in stanene. 

Furthermore, when the size gets smaller than the nanostructuring-induced phonon 
scattering becomes dominant over the three-phonon scattering, and the small-grain-limit 
reduced k becomes proportional to a constant value We calculate the ratio of the 

thermal conductivity to the thermal conductivity per unit of MFP in the small-grain limit; 
for stanene the Isg is found to be 195.0 nm, 92.1 nm, 60.9 nm, 45.6 nm and 36.5 nm at 150 K, 
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FIG. 9: (a) Cumulative lattice thermal conductivity of stanene as a function of the phonon MFP 
at 150 K, 300 K, 450 K, 600 K and 750 K. The curves fitted by Eq. Q are plotted with dot lines, 
(b) Thermal conductivities of stanene nanowires along [100] direction as a function of width. 

300 K, 450 K, 600 K and 750 K, respectively. This quantity is crucial to the thermal design 
for modulating the thermal conductivity in the small-grain limit, for example nanowires. 

It is worthy to mention that for a nanowire system, phonons with long MFPs will be 
strongly scattered by the boundary, thus the contribution from which to k will be limited. 
As shown in Fig. [^b), k decreases with decreasing width of stanene nanowires, and drops 
to half the maximum n in the thermodynamic limit at widths about 1.3 pm. 

In summary, we predict the lattice thermal conductivity of stanene using hrst principle 
calculations and an iterative solution of the BTE. A much lower thermal conductivity is 
observed in stanene than other 2D materials. Our results indicate that stanene, as a 2D 
topological insulator with low thermal conductivity, can realize much higher thermoelectric 
efficiency and is a promising candidate for next-generation thermoelectric devices. The 
predicted contribution of LA phonons to thermal conductivity in stanene is larger than 57% 
due to high phonon group velocities and restricted phase space for three-phonon processes. 
The representative MFP of stanene is obtained for the purpose of future design in TE- 
application nanostructures. The MFP in the small-grain limit is also obtained for designing 
TE devices when the size decreases and the nanostructuring-induced phonon scattering 
dominates. 
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Methods 

All the calculations are performed using the Vienna ab-initio simulation package (VASP) 
based on density functional theory (DFT)E2l. We choose the generalized gradient approx¬ 
imation (GGA) in the Perdew-Burke-Ernzerhof (PBE) parametrization for the exchange- 
correlation functional. We use the projector-augmented-wave potential with Ad electrons of 
tin described as valence, and a plane-wave basis set is employed with kinetic energy cutoff 
of 500 eV. A 15x15x1 k-mesh is used during structural relaxation for the unit cell until 
the energy differences are converged within 10“® eV, with a Hellman-Feynman force conver¬ 
gence threshold of 10“^ eV/A. To eliminate interactions between the adjacent supercells a 
minimum of 15 A vacuum spacing is kept. 

In the calculation of phonon dispersion, the harmonic interatomic force constants (IFGs) 
are obtained using density functional perturbation theory (DFPT), which calculates the 
dynamical matrix through the linear response of electron densitj^. The 5x5x1 supercell 
with 5x5x1 k-mesh is used to ensure the convergence. The phonon dispersion is obtained 
using the Phonopy code with the harmonic IFGs as inpulP^. 

To calculate thermal conductivity of stanene, a n h armonic third order IFGs are needed, 
besides the harmonic IFGs mentioned above. The same 5x5x1 supercell and 5x5x1 k-mesh 
are used to obtain the anharmonic IFGs, and an interaction range of 8.5 A is considered 
herein, which includes hfth-nearest-neighbor atoms. We calculate the thermal conductivity 
of naturally occuring and isotopically pure stanene by solving the phonon Boltzmann trans¬ 
port equation (BTE) using the ShengBTE cod d^^*^^ *^, which is completely parameter-free 
and based only on the information of the chemical structure. 
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